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We discuss the effects of interlayer hopping and the resulting fc^-dispersion in the cuprates within 
the framework of the one-band tight binding (TB) model Hamiltonian. Specific forms of the dis- 
persion relations in terms of the in-plane hopping parameters t, t' , t" and t'" and the effective 
interlayer hopping tz in La2-a;Sra;Cu04 (LSCO) and Nd2-iCe2,Cu04 (NCCO) and the added intra- 
cell hopping tu between the Cu02 bilayers in Bi2Sr2CaCu208 (Bi2212) are presented. The values 
of the 'bare' parameters are obtained via fits with the first principles LDA-based band structures 
in LSCO, NCCO and Bi2212. The corresponding 'dressed' parameter sets which account for cor- 
relation effects beyond the LDA are derived by fitting experimental FS maps and dispersions near 
the Fermi energy in optimally doped and overdoped systems. The interlayer couplings tz and tbi 
are found generally to be a substantial fraction of the in-plane hopping t, although the value of 
tz in NCCO is anomalously small, reflecting absence of apical O atoms in the crystal structure. 
Our results provide some insight into the issues of the determination of doping from experimental 
FS maps in Bi2212, the role of intercell coupling in c-axis transport, and the possible correlations 
between the doping dependencies of the binding energies of the Van Hove singularities (VHSs) and 
various prominent features observed in the angle-resolved photoemission (ARPES) and tunneling 
spectra of the cuprates. 



I. INTRODUCTION 

A good deal of the existing literature on the cuprates 
invokes the one-band tight binding (TB) model Hamil- 
tonian based essentially on the properties of a single 
Cu02-layer. Such a two-dimensional (2D) treatment ig- 
nores intercell coupling in a single layer material such 
as La2-a;Sr:,Cu04 (LSCO) or Nd2_:rCe^Cu04 (NCCO) 
and the additional intracell couplings in multi-layer ma- 
terials like Bi2Sr2CaCu208 (Bi2212). The interlayer hop- 
ping however is responsible for inducing finite dispersion 
of energy bands with kz and for splitting the two Cu02 
bands in Bi2212 into bonding and antibonding combina- 
tions whose existence has become widely accepted in the 
high-Tc community in the last few yearsiiSi^. Moreover, it 
has been shown recently that the angle-resolved photoe- 
mission (ARPES) spectra from a quasi 2D material differ 
fundamentally from the three-dimensional (3D) case due 
to the loss of kz selectivity"*. The ARPES peaks as a re- 
sult display an irreducible width which does not have its 
origin in any scattering mechanism. Such /c^ -dispersion- 
induced linewidths possess characteristic dependencies 
and have been shown to be quite significant in LSCO'^ 
and Bi2212'*, complicating the interpretation of spectro- 
scopic data. 

It is clear with this backdrop that it is important to 
incorporate interlayer couplings and the associated kz- 
dispersion effects in modeling the cuprates. This arti- 
cle represents the first such comprehensive attempt by 
considering the examples of the single layer LSCO and 
NCCO and the double layer Bi2212 systems. For this 
purpose, the framework of the standard dispersion rela- 
tion, E{k) = E[k\\), is generalized to include an addi- 
tional term Ez{k\\ ,kz). Our focus is on the form of Ez in 



LSCO, NCCO and Bi2212 in terms of the effective inter- 
cell, interlayer hopping parameter tz and in the case of 
Bi2212 the added intracell hopping t^i, although forms 
in LSCO and NCCO involving additional intercell hop- 
ping terms are also presented. The values of various TB 
parameters are obtained via fits to the first principles lo- 
cal density approximation (LDA) based band structures 
as well as to the experimental Fermi surfaces (FSs) and 
dispersions near the Fermi energy {Ep) in the optimally 
doped and overdoped systems. In particular, we derive 
estimates of the interlayer hopping parameters in terms of 
the bilayer splitting and the characteristic fcy -dependent 
broadenings of the FS features observed in FS maps ob- 
tained via ARPES experiments. Similarities and differ- 
ences between the in-plane and out-of-plane hopping pa- 
rameters so determined in LSCO, NCCO and Bi2212 are 
delineated. These results are used to gain some insight 
into a number of issues of current interest such as the de- 
termination of doping from the experimental FS maps in 
Bi2212, the role of intercell coupling in c-axis transport, 
and the possible correlations between the doping depen- 
dencies of the binding energies of the Van Hove singular- 
ities (VHSs) and various prominent features observed in 
ARPES and tunneling spectra of the cuprates. 

It should be noted that our analysis implicitly assumes 
the presence of large FSs and an LDA-type electronic 
spectrum, albeit with some renormalization, and it is 
limited therefore generally to the description of the op- 
timally doped and overdoped cuprates. The situation 
with underdoping becomes more complicated, especially 
in LSCO and NCCO, as a lower Hubbard band with 
an increasing spectral weight develops, which is not ac- 
counted for by the LDA or the present TB model fits. 
However, the hopping parameters of the sort obtained 
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FIG. 1: Schematic crystal structure of cuprates, showing Cu 
positions only, to illustrate various hopping processes, (a) In- 
plane hopping parameters t, t' , t" and t'" between different 
Cu atoms in the plane, (b) Simplified out-of plane hopping 
parameters along the c-axis. t-u illustrates hopping within a 
bilayer in Bi2212 and t'^ across the unit cells. The actual 
hopping process is more complicated as discussed in the text. 
Note that the stacking of Cu02 layers across adjacent unit 
cells is staggered, i.e. shifted by (a/2, a/2). 

here based on the "uncorrelated" LDA-type spectrum 
form the traditional starting point for investigating the 
effects of strong correlations and arc thus of consider- 
able value in exploring the physics of the cuprates more 
generally. 

An outline of this article is as follows. Section II 
presents relevant forms of energy band dispersions within 
the framework of the one band model as applied to 
LSCO, NCCO and Bi2212. Section III briefly describes 
the procedure used to obtain fitted TB parameters. Sec- 
tion IV considers the 'bare' parameters which fit the 
LDA-based band structures near the Fermi energy and 
the quahty of the resulting fits. Section V addresses mod- 
ifications needed to fit experimental FSs and dispersions 
and the associated renormalized or 'dressed' parameters. 
Section VI summarizes various bare and dressed param- 
eter sets and comments on other such datasets available 
in the literature. Section VII briefly discusses a few illus- 
trative examples of applications of our results. Finally, 
Section VIII presents a summary and makes a few con- 
cluding remarks. 

II. MODEL DISPERSIONS 

In developing model dispersions for various compounds 
within a one band (single orbital on Cu) TB scheme, we 
begin by recalling that in general the three-dimensional 
energy bands for the TB Hamiltonian can be expressed 
as (assuming a lattice without a basis): 

£:(k) = ^ exp{ik ■ \)ti , (1) 
1 

where 1 denotes a lattice vector and t\ is the hopping 
parameter associated with the 1*^ site. Since our focus 
is on one-band models, we highlight the arrangement of 



the Cu atoms within the Cu02 planes as well as between 
different Cu02 planes in Fig. 1, where various intraplane 
and interplane hoppings are identified. 

In order to emphasize the effect of fc^-dispersion, we 
write 

i?(k) = i?||(k||)-f i?,(k||,fc,). (2) 

Here, k|| and kz respectively denote the in-plane and 
out-of-plane components of k. ^Ey is the model disper- 
sion associated with a single Cu02 plane if the effects 
of fcz-dispersion are totally neglected. The term ac- 
counts for the effect of coupling between different Cu02 
planes, and possesses a very different form in a single 
layer cuprate such as LSCO or NCCO than when a Cu02 
bilayer is present as in Bi2212. We now discuss the two 
terms on the right hand side of Eq. |21 in turn. 

A. In-plane term i5 

i^ll is straightforwardly expressed via the hoppings i, 
t', t" and t'" of Fig. 1(a) as 

E\\ (k||) = -2t[cx{a) + Cy{a)] - At'cx{a)cy{a) 
-2t"[c^(2a) + Cy{2a)] 
-2t"'[c^i2a)cy{a) + Cy{2a)c^ia)] , (3) 

where 

Ci{aa) = cos{akia) , (4) 

and a is an integer (or half-integer) . The shorthand nota- 
tion of Eq.^is convenient as well as transparent since fac- 
tors of cos{akia) occur frequently in describing hoppings 
between lattice sites separated by different distances aa. 
Only a values of 1 and 2 occur in Eq. ^because hopping 
up to a maximum of two lattice sites in the x and/or y 
direction is included as seen from Fig. 1(a). 

B. Inter layer term Ez 

The form of the interlayer component depends sen- 
sitively on the stacking of the layers. Consider an iso- 
lated bilayer first. From Fig. 1(b), the simplest vertical 
hop is straight up via ty . Such a term however in Eq. |21 
will be independent of k\\ . Realistic band computations 
show that this is not the case and that the bilayer split 
bands in Bi2212 display considerable dispersion with fcy. 
In particular, Ref. 6 adduced 

Em = ~2tbi{cj:{a) - Cyia))"^ . (5) 

It should be emphasized that tbi here is an effective pa- 
rameter for hopping within a single bilayer. By expand- 
ing the right side of Eq. 5 as 

{cx{a) - Cy{a))'^ = 1 - 2c3:{a)cy{a) + [ca;(2a) -|- Cj,(2a)]/2 , 

(6) 
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it is easily seen that Ebi in fact involves the sum of three 
different contributions and that its fc|| dependence is the 
result of various in-plane hoppings via the factors of Cx 
and Cy. [The first term of unity represents the vertical 
hop, ty, which is independent of fc||.] Ref. i7' shows how 
form [3 for the bilayer splitting arises in the band struc- 
ture of Bi2212 through hopping from Cu dx2_y2 to Cu 4s 
via the planar oxygen orbitals, with Cu 4s states playing 
a dominant role. 

1. LSCO 

In LSCO, the dominant intercell hopping has a form 
similar to that of Eq. modified by the body-centered 
tetragonal structure, which leads to an offset of successive 
Cu02 planes by half a unit cell in the diagonal in-plane 
direction as seen in Fig.^b), resulting in 

E, = -2t,c,{c/2){cx{a) - Cy{a) fSxy , (7) 

where c denotes the lattice constant along the z-axis, and 

S^y=Cx{a/2)cy{a/2) (8) 

accounts for the presence of the aforementioned layer 
offset^. The term Cz arises because we have an infinite 
number of stacked layers, whereas the expression of Eq.El 
refers to the case of an isolated pair of layers. Note that 
Sxy vanishes for in-plane momenta along (tt, 0) — > (tt, tt), 
leading to a lack of /c^-dispersion along this high symme- 
try line, which is consistent with the observed dispersion 
in LSCO. 



2. NCCO 

In NCCO, the tz term similar to Eq.[Z|is anomalously 
small, due presumably to the absence of the apical oxy- 
gens. Therefore, a number of other competing terms 
come into play. We empirically found that the follow- 
ing form approximates the first-principles band structure 
results: 

Ez = -2cz{cl2)[tz{cx{a) - Cy{a))2 + 2t'^s^X(^l2)sy{a/2) 

+tl{sx{2a) + Sy{2a))]Sxy ,(9) 

where 

Si{aa) = sin{akia) , (10) 
which is similar to the definition of Ci in Eq. 0] 

3. Bi2212 

Bi2212 involves two independent interlayer hopping 
parameters: tbi which controls the intracell bilayer split- 
ting, and tz, which involves intercell hopping. The asso- 
ciated dispersion can be written as 

Ez = -Tz{k\\,Cz{c/2))[{cx[a) - cy{a)f /^ + ao] , (11) 



where 

Tz = ±^Jtl+A'^ + 2tuKcz{c/2) , (12) 

plus (minus) sign refers to the bonding (antibonding) so- 
lution and 

a; EE UzSxy . (13) 

Note that the {cx{a) — Cy{a)Y term leads to absence of kz- 
dispersion along the T — > (tt, tt) line, and particularly at 
r. The introduction of the additional, 'vertical' hopping 
parameter, ao allows for the presence of a splitting at F. 

III. FITTING PROCEDURE 

We first fit the first principles band sructures in vari- 
ous compounds with the TB forms of Eqs.|21-Elto obtain 
what may be referred to as the 'bare' hopping param- 
eters. For this purpose, self-consistent band structure 
calculations were carried out within the LDA using the 
well-established Green's function methodologj^^-. The fit- 
ting process starts by considering the in-plane dispersion 
and can be illustrated with the example of the monolayer 
system. Specifically, the in-plane dispersion is analyzed 
in terms of Eq. |2|by fixing /c^ = tt/c and using the LDA 
energies at a set of high symmetry points. The zero of 
energy is set for convenience at the energy of the VHS 
given by the (tt, 0) point, i.e., 

Ev=E{T:,Q)=A{t'-t"). (14) 

The energies at four other symmetry points from Eq. |31 
then are (with respect to Ey) 

El = E{T) = -A{t + 2t" + t'") , (15) 
E2 ^ E{Tr, tt) = 4(i - 2t" + t'") , (16) 
E3 = E{tt/2,tt/2) ^ A{2t" -t') . (17) 

Ei = E{-K, 7r/2) = 2{t - t'") + 4(i" - t') . (18) 

Fas. 1151 - [THl can be inverted straightforwardly to obtain 
the four unknowns, t,t',t" and t'" in terms of the four 
LDA energies Ei — E4. In terms of the bandwidth 

W = E2- Ei=8{t + t"') (19) 

we have 

t^{Ei + Ev + W/4)/A , (20) 

t' = {W/2-E2)/8, (21) 

t" = t'/2 + E3/8 , (22) 



4 




r (71,0) (71,71) r 




FIG. 2: Comparison of first principles energy bands in LSCO 
(solid lines) with the corresponding tight binding fit based on 
Eq|7|(dashed lines), (a) Dispersion as a function of in- plane 
momentum fc|| along several high symmetry lines for fez = 
and fez = 27r/c. Upper set of bands is shifted by 1 eV with 
respect to the lower set for clarity. E{t:, 0), which is the energy 
of the VHS, defines the energy zero, (b) Effect of fez -dispersion 
as a function of fey for LDA bands. Shaded region is bounded 
by the bands for kz — and fez = 27r/c. (c) Same as (b), 
except that this figure refers to the corresponding results for 
the fitted bands. 
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FIG. 3: FS maps as a function of fey for different fez values in 
LSCO are compared for the LDA (solid lines) and the corre- 
sponding tight binding fit (dashed lines) using Eq. |7| Results 
for doping levels of a; = (a) and x = 0.33 (b) are shown. 



can be modified or 'dressed' to fit the experimental data 
in optimally doped and overdoped systems which exhibit 
large hole-like FSs. The details of the fitting procedure 
vary somewhat for the different compounds considered 
and are discussed below. While the values of the dressed 
and bare parameters differ significantly, the main effect 
appears to be a renormalization (decrease) of the band- 
width by a multiplicative factor Z, with little effect on 
the shape of the FS. 



IV. BARE PARAMETERS 

In this section we discuss the bare parameters and the 
extent to which they reproduce the LDA band structures 
and FSs. 



and 

t'" ^ W/8 - t . (23) 

The resulting dispersions do not completely agree with 
the LDA, indicating that more distant hops are signifi- 
cant. Instead of including these, we have adjusted the 
parameters slightly to optimize the overall fit. In Bi2212 
we apply the above procedure to the average of the two 
bilayer split bands. Once the in-plane parameters are de- 
termined, the out-of-plane {kz) dispersion is fitted along 
similar lines using Eas.lTHTSlto obtain the interlayer hop- 
ping parameters. 

It should be noted that the bottom of the valence band 
in the cuprates usually overlaps with other bands derived 
from planar and apical 0-atoms. This makes it difficult 
to locate the position of the band bottom Ei in LDA 
band structures. For this reason, we have used an average 
value of LDA levels in the vicinity of the band bottom to 
define Ei used in Eq. ^1 

The measured FSs and band dispersions in the 
cuprates in general do not agree with the LDA results, 
indicating the presence of electron correlation effects be- 
yong the LDA. The bare TB parameters discussed so far 



A. LSCO 

The fits in LSCO based on Eqs.Elanddyield: t = 0.43 
eV, t'/t = -0.09, t"/t = 0.07, t"'/t = 0.16 and t^/t = 
0.12. The interlayer hopping tz is a significant fraction 
of the in-plane nearest neighbor (NN) hopping t and the 
size of tz is comparable to the values of in-plane second 
and higher neighbor terms. Fig. [21 considers the quality 
of the TB fit. The LDA and TB bands as a function of fcy 
a,t kz = as well as a.t kz — 2%/ c are seen from Fig. [2Ia) 
to be in good accord near the region of the VHS (defined 
as energy zero) , which is the region of greatest interest in 
delineating near-FS physics. The fc|| dependence of the 
kz dispersion and how it differs between the TB and the 
LDA is better seen by comparing (b) and (c), where the 
shaded regions give the energy range over which the band 
wanders as a function of kz- 

Fig. 131 considers FS maps. Given the good fit to the 
LDA bands near the VHS noted above, the agreement 
between the TB and LDA FSs is to be expected. The 
TB fit correctly reproduces the cuts through the FS at 
kz = and kz = 27r/c as well as changes in the shape of 
the FS as one goes from the undoped case in (a) to the 
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FIG. 4: Same as Fig. I^a), except that the TB fit (dashed 
lines) here is based on Eq. 1241 which better reproduces the 
LDA result (solid lines) near F. 

highly doped case ol x — 0.33 in (b). 

Notably, the TB fit in Fig.[51deviates significantly from 
the LDA bands around F for kz = 27r/c. This is not 
surprising since, as already pointed out, the band bot- 
tom overlaps with other Cu-derived bands. Therefore, 
a proper solution to this problem requires the use of a 
multi-band TB model. Nevertheless, we have empiri- 
cally found that the following form within the single band 
model gives a reasonably good fit. 

Ez = -2(t,c,(l/2)-fi;c2(l/2))[(c,(l)-c,(l))2-Hao5^,]5., 

(24) 

where an additional hopping parameter t'^ has been in- 
troduced together with the constant oq; other quantities 
here have been defined earlier. The fit based on Eqs. El 
and El gives: t = 0.4 eV, t' jt = -0.125, t" /t = 0.05, 
t"'/t = 0.125, tjt = 0.125, t'Jt = 0.05, and = 0.083. 
By comparing these values with those given above based 
on Eqs. El and we see that the in- plane parameters are 
quite similar in the two sets of fits and are thus quite 
robust. Fig.Elshows that the fit based on Ea.l24lcorrects 
the discrepency seen in Fig. El near the F— point. 

Although Eq. provides a better fit to the LDA 
bands, we should keep in mind that there is no clear 
evidence for the presence of an array of bands around a 
binding energy of ~ 1 eV in the cuprates. It is possi- 
ble that electron correlation effects lift the (1^2 _y2 band 
above the complex of Cu- and 0-related bands. If so, 
Eq. El may be better suited than Eq. (31 for describing 
the physical systemii^i 

B. NCCO 

In NCCO, fits were made by using Eq.Elas well as the 
more complicated Eq. El for ■ The resulting in-plane 
parameters are: t = 0.42 eV, t' /t = -0.24, t" /t = 0.15, 
and t'" /t = 0.04. For interlayer parameters, tz/t ~ —0.02 
for both fits, and the additional parameters of Eq.Elare: 
t'Jt = -0.025 and t'Jt = 0.005. The interlayer couphng 
is much smaller than in LSCO, refiecting the absence 
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FIG. 5: (a) Same as Fig. Ela), except that this figure refers 
to NCCO using Eq.Q (b) Effect of fc^-dispersion is shown as 
a function of in the form of the spread E{kz = 0) — E{kz = 
2tt/c) for LDA (solid line) and fits based on Eqs. |7| (dotted 
line) and 1^ (dashed line). 



of apical 0-atoms, as already noted. Figure Eta) shows 
that the TB fit to the LDA bands is quite good. The 
differences between the TB and the LDA bands are at 
the meV level and are displayed in the plot of AEz = 
E{kz = 0) - E{kz = 27r/c) in Fig.E^b). The simple form 
of interlayer hopping in Eq. El (dotted line in Fig. Elb)) 
explains only part of the fc^-dispersion, and the added 
terms in Eq. El (dashed line) improve the fit. The fitted 
FSs based on either Eq. El or Eq- El £^re in very good 
accord with the LDA and are not shown in the interest 
of brevity. 



C. Bi2212 

The fits on Bi2212 are based on Ea.lllland the derived 
parameters are t = 0.36 eV, t'/t = -0.28, t" /t = 0.1, 
t"'/t = 0.06, tu/t = 0.3, tz/t = 0.1, and ao = 0.4. The 
coupling tbi within a bilayer and the intercell coupling 
tz are both quite substantial as are the in-plane hopping 
terms beyond the NN term. Figure El shows that the TB 
fit reasonably reproduces the bilayer split LDA bands de- 
spite some discrepancies. The fits for the FS (not shown) 
are also quite good. 

The lower pair of bilayer split LDA as well as TB 
bands in Fig. El for kz = show an interesting level 
crossing between the B (bonding) and A (antibonding) 
bands around k = k| = (0, 0.47r/a, 0). This level crossing 
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FIG. 6: Same as Fig. |2Ia), except that this figure refers to 
Bi2212 and is fit using Eq. 1111 Note that the bilayer spUt 
bonding (B) and anti-bonding (A) bands are shown for each 
fcz value. 



is associated with an anomalous /cx-dispersion, which is 
highlighted in Fig.[71 the TB bands are drawn on an ex- 
panded scale along the F— (tt, 0) line in (a), while the k^- 
dispersion at three illustrative values is shown in (b)- 
(d). The shape changes from being described approxi- 
mately by the form sin^(A:2c/4) for fcy = 0, to | sin(A:2c/4)| 
at kZ, and finally back to sin^(fczc/4) for fc||a/(7r) = 1.0. 




(7t,0) 



> 

o 

IW 

W 



0.0 



-0.5 



1 ■ 1 ■ 1 ■ 

'(b) 


'(cj ' ■ ' ■ 








X=0-437i;/a" 









k^=0.907i:/a 



J I I I L 



0.0 0.8 0.0 0.8 0.0 0.8 

FIG. 7: (a) Effect of fc^-dispersion for the TB fit in Bi2212 
showing the crossing between the bonding and antibonding 
bands around fcx = 0.437r/a. Shaded areas are bounded by 
bands for k\\ along the V — X line for fcz = (dashed lines) and 
fcz = 27r/c (solid lines), (b)-(d) Dispersion of bands with kz 
at three different values of k^^ (shown as vertical dotted lines 
in (a)) showing change in the form of fcz-dispersion through 
the band crossing point in (a). 



V. DRESSED PARAMETERS 

In this section we discuss modified or 'dressed' TB pa- 
rameters needed for describing the measured FSs and 
dispersions. As already pointed out in the Introduction 
above, our analysis here implicitly assumes the presence 
of large LDA-type FSs, and it is thus limited generally 
to the optimally doped and overdoped systems. 

We emphasize that comparisons between the measured 
and computed FS maps only yield information on the 
relative sizes of the hopping parameters, e.g. the ratios 
ti/t for various U^s in terms of the NN hopping t, and 
that dispersions are additionally needed to derive the en- 
ergy scale t. Unfortunately, experimental ARPES data 
on dispersions are typically available only over a limited 
binding energy range of about 1 eV, leading to inherent 
uncertainty in any fitting process. Our procedure for ob- 
taining the dressed parameters starts by comparing the 
experimental FSs with the ones computed by using the 
bare parameters of the previous section. In this connec- 
tion, we compute FS imprints in the (fc^;, fcj,)-plane for 
kz — and kz = Stt/c, in order to identify the momen- 
tum region over which the ARPES intensity is allowed 
due to the effect of fcz-dispersion along the lines of Refs.0 
and^;. In this way, the interlayer coupling parameters can 
be estimated. All parameters are then adjusted to obtain 
a good fit. Fig. 8 gives insight into these results and the 
quality of the final FS fits in the three systems considered. 
Finally, the measured dispersions are fitted to determine 
the average value of the renormalization factor Z. Fur- 
ther relevant specifics in the individual cases of LSCO, 
NCCO and Bi2212 are given in the Subsections V A-C 
below. In the interest of brevity, the final values of all 
dressed parameters are listed in Table II below. 



A. LSCO 

In LSCO, the experimental FS maps are in good over- 
all accord with the first principles ARPES computations^ 
over an extended doping range of a; = 0.06 — 0.22, where 
in the underdoped regime, this agreement suggests that a 
remnant of the LDA-type quasiparticles with a reduced 
spectral weight continues to persist in the presence of 
strong correlations. In this work, we have specifically 
chosen to fine tune the bare TB parameters to fit the mea- 
sured FS mapiiiS a.t X = 0.15 using Eq.[3 Fig. i^a) is 
relevant in this regard. The left hand side shows the mea- 
sured ARPES intensity from the FS for reference. The 
right hand side shows the same intensity, except that we 
have overlaid the FS cross sections in the (fca;, fcj,)-plane 
for kz = (solid thick lines) and kz — 27r/c (dashed 
thick lines) obtained by using the final values of the 
dressed parameters (see Table II below). The parame- 
ter ratios are the same as in the LDA fit, except for a 
larger t'" /t — 0.21. The experimental intensity is seen to 
be reasonably well contained within the boundaries de- 
fined by these FS imprints given by the thick solid and 
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FIG. 8: Left hand column gives experimental FS maps (i.e. 
ARPES intensities for emission from the Fermi energy) for: 
(a) An optimally doped LSCO sample with x — 0.15—*—; (b) 
An optimally doped NCCO sample with x = O.lSi^iii; (c) An 
overdoped Bi2212 sample— (identified as OD70 in Ref. Uyt) . 
Right hand side panels give the same spectra together with 
the computed cross sections of the FS in the (fc^,, fcy)-plane 
for kz — (thick solid lines) and kz = 27r/c (thick dashed 
lines); see text for the details of the dressed tight-binding 
parameters used. Note, the solid and dashed FS contours are 
well separated for the LSCO and the antibonding FS sheet 
(A) in Bi2212, but this is not the case for NCCO and the 
bonding sheet (B) of the FS in Bi2212. 

dashed lines. [The ARPES intensity in this region of al- 
lowed transitions is not expected to be uniform as it gets 
modulated considerably in general by the ARPES matrix 
element.] In particular, the value of the dressed parame- 
ter tz used reproduces the characteristic shape of the ex- 
perimental FS map. Our estimate of the tz so obtained 
should be considered an upper bound since the effects of 



experimental resolution and various possible scattering 
mechanisms not accounted for in our model will serve to 
further broaden the theoretical FS imprint. Finally, we 
have fitted the measured dispersioni^, to adduce an av- 
erage renormalization factor of Z = 0.64. The individual 
renormalization factors are Z{t) = Z{t') = Z{t") = Z{tz) 
=0.6, and Z{t"') = 0.8. 



B. NCCO 

In NCCO, the ARPES experimentsi^iSJML seem to 
be dominated by the presence of upper and lower Hub- 
bard bands, which gradually collapse with doping^^. 
Near optimal doping the splitting between the upper and 
lower Hubbard bands is small and the resulting FS map 
resembles the LDA results. Accordingly, we have used 
the X — 0.15 dataset shown in Fig. 8(b) in the fitting 
process based on Eq. 7. However, the size of the mea- 
sured FS is best described if the doping is taken in the 
computations to be a; = 0.16, and therefore, this slightly 
modified value of doping is used in the present fits. Also, 
the observed FS curvature is better fitted by reversing 
the sign of the bare t'" (see Tables I and II below). Its 
effect is to move the FS closer to the magnetic zone 
boundary given by the zone diagonal, and thus this sign 
change may be related to the presence of some residual 
anti-ferromagnetic order in the system. Aside from this 
change, the fitted in-plane parameter ratios are the same 
as the corresponding LDA ratios. 

In discussing the bare parameters in Section IV B 
above, we have already pointed out that the interlayer 
coupling tz in NCCO is quite small. As a result, changes 
in the cross section of the FS with kz are also quite small 
and for this reason the two FS imprints are not resolved 
on the right hand side of Fig. 8(b) (thick solid line). The 
situation is in sharp contrast with the case of LSCO in 
Fig. 8(a) where the kz ~ and kz — 27r/c FS cross- 
sections are well separated. These remarks will make it 
clear that it is difficult to determine the dressed value of 
tz in NCCO with any accuracy via comparisons of the 
measured and computed FSs, and therefore, no value for 
tz is given in Table II below. 

The TB fit to the FS is seen from Fig.EIb) to be quite 
good overall. The experimentally observed linewidths in 
ARPES spectra in NCCO are generally larger than the 
linewidths that can be associated with the effect of kz- 
dispersion. This suggests that various scattering mecha- 
nisms are at play especially around the 'hot spots' where 
the LDA FS crosses the zone diagonal, allowing for the 
opening of pseudogaps which are not accounted for in 
the LDA or the present TB fits. As in the case of LSCO, 
the final values of the dressed parameters in NCCO also 
are more or less uniformly renormalized with respect to 
the bare values with a renormalization Z = 0.55 for all 
parameters, except t'" as noted above. 

The present TB parameters can be compared to values 
extracted from strong coupling calculations. In NCCO, 
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FIG. 9: Experimental FS maps obtained via ARPES mea- 
surements by Ref.'22'for six different samples of Bi22I2 cover- 
ing an extended doping range from underdoped (UD) to over- 
doped (OD) are shown overlaid with the corresponding TB FS 
maps; see text for details of the TB parameters and the Fermi 
energies used for different doping levels. Computed bonding 
and antibonding FSs are shown only in the right portion of 
each of the six FS maps by black lines where the thickness of 
the lines represents the width of the FS induced by the effect 
of fcz- dispersion. 

the ARPES datski^ were fit to an antiferromagnetic Slater 
modeU^i. The hopping parameters found from the fit 
are comparable to the bare parameters with an average 
renormalization Z = 0.7 {Z{t) = 0.77, Z(t') = 0.84, 
Z{t") = 0.43). 



C. Bi2212 

The right hand side of Figure (Sfc) shows the TB fit 
based on Eq. ^2 to the experimental normal state FS 
map of an overdoped Tc = 70 K Bi2212 sample, for which 
the bilayer splitting is clearly resolvedi^ii^. For the an- 
tibonding FS sheet, the experimental intensity is seen to 
fall quite well within the fcz = and fc^ = 27r/c slices 
of the computed FS, especially near the antinodal re- 
gion. On the other hand, the bonding FS sheet displays 
little width associated with fc^-dispersion, which is also 
consistent with the character of the measured spectra. 
The fit in Fig. |SIc) was used to determine the dressed 
parameter ratios, while the overall renormalization fac- 
tor was found by fitting the dispersion near (tt, 0)^2*21. 
The renormalization factors for various parameters are 
found to be: Z{t) = 0.35, Z{t') = 0.36, Z{t") = 0.43, 
Zit'") = 0.15, Z{tu) = 0.27, and Z(i^) = 0.14. The 
renormalization of the hopping parameters is thus seen 
to be quite substantial with an average value of Z = 0.28. 
t\,i which is responsible for the bilayer splitting is reduced 
to about 27% of its bare value. Notably, the experimen- 
tal FS is 'pinched' near (7r/2,7r/2), with the nodal point 
lying closer to the F point than in the LDA calculation. 

Using our dressed TB parameters for Bi2212, we have 
attempted to model the doping dependence of the FS 
of Bi2212 by adjusting the Fermi energy as a function 



TABLE I: Bare TB parameters obtained by fitting the first 
principles LDA-based band sructures as discussed in Sec- 
tion IV. 



Parameter 


LSCO 


NCCO 


Bi2212 


t 


430 meV 


420 meV 


360 meV 


t' 


-40 


-100 


-100 


t" 


30 


65 


35 


t'" 


70 


15 


20 




50 


-8 


36 








110 



TABLE II: Dressed TB parameters obtained from fits with 
the experimental FSs and band dispersions as discussed in 
Section V. Value of dressed in NCCO is undetermined, as 
discussed in Section V B. 



Parameter 


LSCO 


NCCO 


Bi2212 


t 


250 meV 


230 meV 


126 meV 


t' 


-25 


-55 


-36 


t" 


20 


35 


15 


t'" 


55 


-30 


3 




30 




5 


tbi 






30 


z 


0.64 


0.55 


0.28 



of doping. [See Section VII A below for the specifics of 
the Fermi energy values used.] The results are shown in 
Fig. ini where the TB FS maps (broadened to reflect the 
effect of /c^-dispersion) are superposed on the experimen- 
tal data of Ref. 122 for a series of dopings from under- 
to the overdoped samples. Since the data were taken at 
room temperature, the bilayer splitting is not clearly re- 
solved in the experimental plots. However, the TB fits 
are seen to reproduce the overall changes in the shape and 
size of the FS in Bi2212 with doping reasonably well. 



VI. SUMMARY OF TIGHT-BINDING 
PARAMETERS 

We collect various TB parameter sets in this section 
for ease of reference and intercomparisons between the 
resuhs for LSCO, NCCO and Bi2212. Table I lists bare 
values, while the dressed values are given in Table II. The 
LSCO datasets in Tables I and II refer to fits based on 
Eq. El although a fit using Eq. 1241 was also discussed in 
Section IV A above. Table III lists a few other parameter 
sets in these compounds available in the literature. 

Some comments concerning the bare parameters in 
LSCO, NCCO and Bi2212 as well the corresponding 
dressed values have already been made in Sections IV 
and V above and need not be repeated. Here we make 
a few further observations as follows. In the single layer 
systems LSCO and NCCO the dressed values are in ex- 
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TABLE III: Other TB parameter sets determined from LDA band structure calculations (columns 1 and 2) or from ARPES 
experiments (columns 3-6). 





LSCO 


Bi2212 


general 


Par am. 


Ref. 8 


Ref. 8 


Ref. 26 


Ref. 27, 


Ref. 28 


Ref. 2^ 


t 


420 meV 


420 meV 


150 meV 


150 meV 


390 meV 


350 meV 


t' 


-73 


-110 


-40 


-24 


-82 


-120 


t" 


36 


55 


13 


33 


39 


80 


t'" 






28 


13 












-13 


-24 
















90 





cellent agreement with the bare values, differing only by 
an overall renormalization Z except for the most distant 
neighbor t'" . For Bi2212 bare and dressed parameters 
have larger differences (associated with the bilayer split- 
ting) and the average renormalization factor is smaller, 
suggestive of stronger correlation effects. Given this gen- 
eral similarity, we restrict our further remarks mainly to 
the bare parameter values. 

Comparing monolayer LSCO and NCCO, the bare 
values of the in-plane NN hopping t are similar (~ 
0.4 eV), but other in-plane hoppings, t' and t" are 
about twice as large in NCCO compared to LSCO. 
For hole-doped cuprates, an increasing value of t' has 
been correlated with increasing superconducting tran- 
sition temperatur o^'^ - ? , and our results for LSCO and 
Bi22f2 are consistent with that trend. The larger t' in 
NCCO suggests that the trend does not hold for electron- 
doped cuprates, but that might be expected since Ep in 
electron-doped cuprates is very far from the VHS. No 
clear trends emerge for t" and t'", except that they are 
comparable in magnitude to t' , and hence should be re- 
tained in detailed analyses. 

One may estimate the strength of strong coupling cor- 
rections from the renormalization factors Z via the rela- 
tion, Z = I/(IH-A), where A is a measure of the com- 
bined strength of electron-electron and electron-phonon 
coupling. This yields A values of 0.6, 0.8 and 2.6 in 
LSCO, NCCO and Bi22I2, respectively, suggesting much 
stronger correlations in Bi22I2, which has the highest 
Tc. These values of A should be used with caution, since 
ARPES tends to underestimate dispersion, even in sim- 
ple metals--, and a more careful evaluation is needed to 
properly account for 'kink' physics^. 

Turning to the nature of interlayer coupling, the values 
of the bare and dressed and ti,i are plotted in FigurelTUI 
as a function of the separation Cz between Cu02 planes. 
It is interesting that this coupling shows a systematic 
decrease with increasing separation between the Cu02 
layers in Bi2212 and LSCO, but that NCCO which has a 
different relative arrangement of Cu and O atoms in dif- 
ferent layers falls outside this pattern. The effect is more 
pronounced for the bare parameters (circles) compared 
to the dressed ones (triangles). 

With regard to other parameter sets for the cuprates 




5 10 

c,(A) 

FIG. 10: Magnitude of the interlayer coupling as a function 
of the layer distance Cz as measured by the values of bare 
(circles) and dressed (triangles) \tz\ and \tbi\ in LSCO, NCCO 
and Bi2212. Lines drawn on the data are guides to the eye. 



in the literature, note that Table III is essentially lim- 
ited to in-plane couplings since little attention has been 
paid to the importance of interlayer coupling in the past. 
The data of Table III are generally similar to the results 
of Tables I and II depending on the type and details of 
fit invoked by various authors. More specifically, Ref. 
(first two columns of Table III) adduces t ~ 0.42 eV for 
a large variety of cuprates, and that t" ^ — <'/2; the for- 
mer result is in agreement with Table I for LSCO, while a 
somewhat smaller value is found for Bi2212. For Bi22I2, 
Refs. |2^ and |23 provide two parameter sets, which in- 
volve more distant neighbors than we have included, al- 
though no values for bilayer or other kz related disper- 
sion are included. [The term reflects a contribution 
— 2i™Ca;(2)cj,(2) to in-plane dispersion.] These parame- 
ter values in Table III are seen to be quite close to our 
dressed values in Table II for Bi22I2. The parameters 
found by Ref. are strikingly different, although the 
parameter ratios are more similar. It may be that these 
parameters refer to the total dispersion, while our evalu- 
ation (Table II) is limited to the region close to the FS. 
It is known that Z is larger in the near-FS regime^. The 
last column of Table III, marked general, lists the param- 
eters reported in Ref. 29 as providing a generalized fit to 
the experimental ARPES data on a variety of cuprates 
and it is difficult to assess their accuracy in the context 
of the individual compounds considered in this work. 



10 



TABLE IV: Estimated values of superconducting transition 
temperature T^, Fermi energy Ep and doping x, for vari- 
ous overdoped (OD), optimally doped (OP), and underdoped 
(UD) Bi2212 samples. 



Rof. No. 


Doping 






X 


19 


OD 


70 K 


-160 meV 


0.29 


21 


OD 


55 


-160 


0.29 


21 


OP 


91 


-130 


0.17 


21 


UD 


78 


-120 


0.135 



VII. ILLUSTRATIVE APPLICATIONS 
A. Doping Dependence of Tc in Bi2212 

The doping level in a given single crystal of Bi2212 
can in principle be determined by measuring the FS and 
the associated volume and hence the hole count via an 
ARPES experiment. ^'^ In practice however the imprint 
of the FS in ARPES becomes broadened due to k^- 
dispersion and other effects, adding uncertainty to such a 
determination of the FS volume in terms of models where 
the interlayer coupling is essentially ignored as is the case 
in much of the existing literature. Here we attempt to 
address this issue anew by using our more complete mod- 
eling of the FS in which the effects of fc^-dispersion and 
bilayer splitting in Bi2212 are explicitly included. 

Our procedure for analyzing the ARPES data in any 
specific Bi2212 sample is as follows. Using the dressed 
TB parameters listed in Table II, we compute the ex- 
pected FS imprints for a range of values of Ep^ including 
the effect of broadening of the FS due to fcz-dispersion. 
[Note that Ep is the only free parameter in the compu- 
tation.] The theoretical FS imprint is then overlaid on 
the experimental one to obtain the value of Ep which 
provides the best fit to the experimental data. Once Ep 
is fixed in this way, the FS volume and the hole doping 
is uniquely determined. 

Using the preceding approach, we have analyzed the 
experimental datasets of Refs. and |^ The results 
of the doping values x obtained for different samples are 
summarized in Table IV. The best fit to i?F is also listed 
in each case. Note that the data cover a substantial dop- 
ing range and involve underdoped (UD), optimally doped 
(OP) as well as overdoped (OD) samples. The Tc vs x 
values of Table IV are plotted in Fig. ^2 and fall rather 
nicely on a parabola. 

In their room temperature measurements, Kordyuk, 
et ali2^ were unable to resolve the bilayer splitting. The 
doping they estimated from the area of the single FS 
was too small, leading them to hypothesize that they 
could account for the effect of bilayer splitting by rigidly 
shifting the data by a constant Ax = 0.07. In Fig. ^2 we 
show that their data when shifted by Aa; = 0.09 line up 
reasonably well with other results. These shifted values 
of X in turn give the theoretical FSs plotted in Fig. |51 




FIG. 11: Tc vs doping x in Bi2212 as determined by using the 
Fermi energies given in Table IV and discussed in the text. 
Experimental data sets are from: Ref . ITgl (square'). I21I (solid 
dots), and Ref. l2^ (diamonds). Also shown is the Ta{x) dome 
for LSCO (dotted line), after Ref.lsJ 

above as broadened black lines. The residual differences 
may indicate subtle differences in the samples^ or that 
the actual doping is not exactly described by a rigid shift. 
These differences also show up in the theoretical FS maps 
in Figs.lHIc) andlHI 

The solid line in Fig. II II gives the parabolic fit through 
all the data points and can be expressed as 

r, = T,o[l-(^^^)^], (25) 

Xl 

with Tco = 95 K being the maximum Tc, a;o = 0.2 is 
the optimal doping, and xi = 0.15 is the width of the 
parabola. This result should allow the conversion of 
data to actual doping values in Bi2212 more generally 
(as long as it is known whether the sample is over- or 
underdoped) . 

Figure [TTl also compares the superconducting dome of 
Bi2212, with that of LSCO^i, which is given by Ea.l^ 
with xq — 0.16, Xl = 0.11, and T^o — 38 K. Both formu- 
las give an onset of superconductivity at a; = 0.05, while 
the doping at optimal Tc in Bi2212 of xq — 0.2 is some- 
what higher. Honma, et al.'^^ have recently reported an 
even higher optimal doping for Bi2212, x ~ 0.25. The 
similarity with the Uemura™ plot is striking, and may 
give insight into why the Uemura plot peaks at different 
dopings for different cuprates^. 



B. c-Axis Conductivity and Intercell Coupling 

Whether the c-axis conductivity CTc in the cuprates is 
coherent or not is currently a subject of considerable 
controversy. However, for a strongly correlated system, 
one would expect that increasing the interlayer hopping 
should increase CTc, whether this conductivity is coherent 
or not. For instance, in a fluctuating stripe model, inter- 
layer conductivity would be sensitive to whether charged 
stripes in different layers are in registry, but the actual 
hopping will be controlled by t^. Also, in the interlayer 
pairing model, single particle interlayer hopping is frus- 
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FIG. 12: c-axis conductivity tJc measured just above Tc 
for optimally doped cuprates plotted vs normalized interlayer 
hopping (tz/t)'^. Filled circles give data of Ref. l3^ solid line 
represents expected form <Jo{tz/t)^ , with ao = 680(f2 — cm)~^ . 




FIG. 13: Solid lines give TB densities of states (DOSs) 
for LSCO and Bi2212, using dressed parameters of Table II. 
Dashed lines give DOSs for the limiting 2D case of tz = 0. 



trated by strong correlation effects, but the imfrustrated 
pair hopping still depends on the bare tz- 

Our purpose here is only to ascertain whether our val- 
ues of tz are consistent vi^ith the observed material varia- 
tions of the c-axis conductivity in the cuprates. In partic- 
ular, Ref. 35 argues that <Tc at fixed doping should scale 
approximately as (tz/t)^. Figure IT^ shows that the plot 
of experimental^^ (Tc vs {tz/t)^ is nearly linear'^^, indicat- 
ing that our tz values are at least not obviously inconsis- 
tent with experiment in LSCO, NCCO and Bi2212. Note 
that we have taken the tXc values in Fig. El for the op- 
timally doped samples where the c-axis conductivity is 
likely to be more coherent. 



C. Effect of kz Dispersion on Van Hove 
Singularities 

Van Hove singularities in the electronic density of 
states (DOS) have been implicated for playing an im- 
portant role in determining various properties of the 
cuprates. If the interlayer coupling is neglected, i.e. the 
system is assumed to be strictly 2D, the VHS occurs as 
a sharp feature at a single energy in the DOS of a mono- 
layer system and at two different energies (correspond- 
ing to the bilayer split bands) in a double layer system. 
Fig. 1131 shows this effect with the examples of LSCO and 
Bi2212. Here the horizontal axis is given in terms of 
the doping x, which can be obtained straightforwardly 
for any given value of the Fermi energy by counting the 
number of holes in the band, allowing the energy scale in 
the DOS to be converted to the doping scale. The single 
sharp VHS peak for LSCO in Fig. [T^a) at x = 0.2 in 
the 2D case (dashed line) is replaced by a fiat region of 
high DOS spread over x = 0.09 — 0.21 when the effect of 
/c^-dispersion is included (solid line). A similar effect is 
seen in Fig. Elb) for Bi2212, except that here we have 
two peaks related to the bonding and antibonding bilayer 
bands. The antibonding VHS (A, lower peak around 
X = 0.32) displays a smaller width compared to the bond- 
ing VHS (B, around x — 0.62), reflecting differences in 
the extent to which A and B bands are spread in en- 
ergy through the coupling tz- By comparing Figs. IH^ b) 
and llll it can be seen that in Bi2212 the optimal Tc cor- 
responds to a doping close to but less that the doping at 
which the antibonding VHS crosses the Fermi level. 

The doping dependence of the positions of the VHSs 
and their possible correlation with superconducting prop- 
erties is briefly explored in Fig. ^1 In order to help ori- 
ent the reader, we consider Fig. I14f a^ for LSCO first. 
The thick solid curve (upper boundary of the shaded re- 
gion) gives the position of the leading edge of the VHS, 
which is seen from Fig. I13r a) to lie at the Fermi en- 
ergy for X = 0.19, and to move to higher binding en- 
ergies with decreasing doping. The dashed thick curve 
similarly gives the binding energy of the trailing edge of 
the VHS. Also shown are the gap parameter A, defined 
via A — 2.05fcsTt,, for the superconducting dome (dot- 
ted line) and the value of the pseudogap from tunneling 
experiments^ (triangles). In the same vein. Fig. I14f b) 
which refers to Bi2212, gives variations in the binding 
energies of the leading and trailing edges of the VHSs for 
the bonding and antibonding bands, together with the 
gap parameter for the superconducting dome and vari- 
ous features in the tunneling and ARPES spectra^i^Si^ 
as detailed in the figure caption. 

Fig. ll3r a) shows that the effect of fcz-dispersion is quite 
substantial on LSCO in that it spreads the VHS over a 
large portion of the superconducting dome. In Bi2212 
on the other hand the effect is smaller and the bonding 
and antibonding VHSs remain distinct. One of the ex- 
perimentally observed features in ARPES and tunneling 
(the lower binding energy datasets in Fig.^Jb)) seems to 
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FIG. 14: Doping dependence of VHS in LSCO (top) 
and Bi2212 (bottom), assuming rigid band filling and using 
dressed parameters of Table II. The shaded features bounded 
by solid lines (leading edge) and dashed lines (trailing edge) 
represent the width of the VHS peaks. Shown for comparison 
is the doping dependence of Tc, plotted as an effective gap, 
A = 2.05kBTc (dotted lines). Various symbols are taken from 
experiment: In (a) (LSCO) downward triangles represent tun- 
neling pseudogaps, Ref.l38t in (b) (Bi2212) tunneling pseudo- 
gaps are represented by upward triangles (Ref. js^; ARPES 
leading edge pseud ogap s by downward triangles (Ref.'39^ and 
filled squares (Ref. l40l) : the ARPES peak by filled diamonds 
(Ref.0); and the ARPES hump by open diamonds CRef. liH) 
and open squares (Ref. |^). 



track the antibonding VHS in Bi2212 over a fairly large 
doping range, but this correlation for the higher bind- 
ing energy feature is hmited to optimal and overdoped 
regimes. In LSCO, the tunnehng feature shown corre- 
lates with the leading edge of the VHS and is well sep- 
arated from the behavior of the trailing edge. However, 
the connection between the tunneling pseudogap (which 
is approximately symmetrical about Ep) and the VHS 
(which is not) is nontrivial, and a more detailed analysis 
is considered beyond the scope of this paper. 



VIII. SUMMARY AND CONCLUSIONS 

We discuss tight-binding parameterization of electronic 
states near the Fermi energy in LSCO, NCCO and Bi2212 
where the effect of fcz-dispersion is included in a compre- 
hensive manner for the first time. Our analysis proceeds 
within the framework of the one-band model and incorpo- 
rates available ARPES data and the results of first princi- 
ples band structure computations on these cuprates. The 



in-plane coupling is treated in terms of the hoppings t, 
t', t" , and t'" between the nearest and farther out neigh- 
bors, while the coupling between different CUO2 layers is 
considered mainly through an effective intercell hopping 
tz , with an added hopping thi to account for the intracell 
interaction involving the Cu02 bilayer in Bi2212. In each 
compound we obtain the specific tight-binding form for 
describing the dispersion and the associated 'bare' pa- 
rameter sets, which fit the 3D band structures based on 
the LDA, as well as 'dressed' parameters which account 
for correlation effects beyond the LDA on the experimen- 
tal FSs and near-FS dispersions in optimally doped and 
overdoped regimes. 

In LSCO, the bare interlayer hopping — 0.05 eV 
is found to be a substantial fraction of the in-plane NN 
hopping t with tz/t = 0.12, and the size of tz is compa- 
rable to the values of in-plane second and higher neigh- 
bor terms. In NCCO, the interlayer coupling is anoma- 
lously small with tz/t ~ —0.02, reflecting lack of apical 
0-atoms in the crystal structure; we delineate how in this 
case other smaller terms become important and modify 
the kz dependence of the dispersion. In Bi2212, the bi- 
layer coupling tbi and the interlayer coupling tz are both 
quite substantial, with tz/t = 0.1 and tbi/t = 0.3, as are 
the in-plane terms beyond the NN term; a level degener- 
acy between the bonding and antibonding bands around 
k = (0, 0.47r/a, 0) induces an anomaly in kz dispersion. 
Correlation effects beyond the LDA reduce the band dis- 
persions as expected, so that values of the dressed pa- 
rameters are smaller than the corresponding bare values. 
Renormalization factors for the hopping parameters vary 
substantially with average values of Z = 0.64, 0.55, and 
0.28 in LSCO, NCCO and Bi2212, respectively. 

It should be noted that any fitting procedure involving 
a limited basis set - a one band model here - possesses 
inherent uncertainty tied to the binding energy used in 
obtaining the fit. In the cuprates, the present frame- 
work of the one band model would need to be modified 
when one moves to higher binding energies and bands de- 
rived from other Cu d- and O p -orbitals come into play. 
Also, the experimental values of the interlayer couplings 
tz and tbi adduced in this work should be viewed as up- 
per bounds since these are based on model computations 
which only include the broadening effect of /c^-dispersion 
and do not take into account the effects of various scat- 
tering mechanisms and experimental resolution. 

Our modeling of the /c^-dispersion effects provides in- 
sight into a number of issues of current interest and we 
briefly illustrate this aspect with three examples, (i) 
We show that the inclusion of fc^-dispersion in model- 
ing FS imprints yields doping levels x for various pub- 
lished ARPES datasets in Bi2212, which when plotted 
against related T^s, fall remarkably well on a parabola. 
This parametrization of the superconducting dome may 
help the determination of doping levels of Bi2212 single 
crystals via ARPES experiments more generally, (ii) By 
using the values derived in this work, we show that the 
intercell hopping controlled by t^ is consistent with c-axis 
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transport, at least in optimally doped samples, with the 
c-axis conductivity scaling approximately as (t^/t)'^- 
(iii) Logarithmic Van Hove singularities (VHSs) in the 
2D case are spread over energy when /cx-dispersion is in- 
cluded. We compare the binding energies of the lead- 
ing and trailing edges of the VHSs with those of various 
prominent features observed in ARPES and tunneling ex- 
periments in LSCO and Bi2212. The antibonding VHS in 
Bi2212 tracks some experimental features quite well over 
a wide doping range, but this is the case for the bond- 
ing VHS in Bi2212 and the leading edge of the VHS in 
LSCO mostly in the optimally doped regime. Although 
specific implications of these results require further anal- 
ysis, they suggest that attention should be given to the 
kz induced width of the VHSs when making comparisons 
with experimental features. 



We conclude that the simple tight-binding 
parametrization of LSCO, NCCO and Bi2212 pre- 
sented in this article constitute a useful starting point 
for exploring the effects of fc^-dispersion on the electronic 
properties of the cuprates. 
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